Transformation rule for covariance matrices under Bell-like detections 
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Starting from the transformation rule of a covariance matrix under homodyne detections, we can 
easily derive a formula for the transformation of a covariance matrix of (n + 2) bosonic modes under 
Bell-like detections, where the last two modes are combined in an arbitrary beam splitter (i.e., with 
arbitrary transmissivity) and then homodyned. This formula can be specialized to describe the 
standard Bell detection and the heterodyne measurement, which are exploited in many contexts, 
including protocols of quantum teleportation, entanglement swapping and quantum cryptography. 
Our general formula can be adopted to study these protocols in the presence of experimental im- 
perfections or asymmetric setups, e.g., deriving from the use of unbalanced beam splitters. 
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I. INTRODUCTION 

Gaussian quantum information (GQI) is that area of 
quantum information which deals with continuous vari- 
able systems (e.g., bosonic systems) prepared in Gaus- 
sian states, evolving by Gaussian unitaries or channels 
and measured by Gaussian detectors [l[ . Gaussian states 
are easy to generate experimentally and easy to manip- 
ulate theoretically. Their description can be reduced to 
their first two statistical moments, which are the mean 
value (or displacement vector) and the covariance matrix 
(CM). In particular, the CM contains the most relevant 
information about the Gaussian state, providing its en- 
tropy, purity properties and separability properties [H, [1J . 

One of the most important Gaussian measurements in 
GQI is the Bell detection [l]. This consists of combining 
two bosonic modes into a balanced-beam splitter (i.e., 
with transmissivity 1/2). The output modes are then 
measured by two homodyne detectors in such a way that 
one mode is detected in the position quadrature q and 
the other mode in the momentum quadrature p. This 
measurement is common in a series of protocols, e.g., 
in continuous variable quantum teleportation [3j-|9| and 
cntaglement swapping [lCHl2j . Another important mea- 
surement is heterodyne detection, where a single bosonic 
mode is taken as input of a balanced-beam splitter (with 
the other input being the vacuum) and the two outputs 
are homodyned in q and p, respectively. This is a fun- 
damental detection in many protocols, for instance in 
coherent-state quantum key distribution [l3l - fl(| and two- 
way quantum cryptography [17| . 

In this paper, we consider a generalized measurement 
that we call Bell-like detection. Here we have two bosonic 
modes which are combined into a beam splitter of ar- 
bitrary transmissivity and then homodyned in the two 
quadratures (one mode in q and the other in p). Stan- 
dard Bell detection and heterodyne detection are spe- 
cific instances of this more general measurement. In our 
derivation, we consider a set of n + 2 bosonic modes in a 
Gaussian state with a generic CM. By applying the Bell- 
like detection to the last two modes, we compute the 
conditional reduced CM of the first n modes surviving 



the measurement. Our algebraic derivation is relatively 
easy starting from the well-known transformation rules 
of CMs under partial homodyne detections [H, [l9| . De- 
spite its easy derivation, our main formula can be usefully 
applied to several contexts. For instance, it can be ex- 
ploited to extend protocols of teleportation and entangle- 
ment swapping to the case of unbalanced beam splitters 
(i.e., imperfect Bell measurements). Similarly, it can be 
used to perturb the ideal heterodyne detection which is 
used in many protocols of quantum key distribution. 

The paper is organized as follows. In Sec. [II] we provide 
a brief introduction to bosonic Gaussian states and CMs. 
In Sec. IIIII we review the transformation rule of CMs un- 
der homodyne detections. Then, in Sec. IIV1 we derive 
the main result of the paper, i.e., the formula for the 
transformation of CMs under Bell-like detections. Spe- 
cific formulas for the standard Bell detection and het- 
erodyne detection are also shown. Finally, Sec. [V] is for 
conclusions. 



II. BASIC NOTIONS ON BOSONIC GAUSSIAN 
STATES 

A system of n bosonic modes is a described by a vector 
of In quadrature operators 



(91, Pi, ■ • -,Qn,Pn) , 



(1) 



satisfying the commutation relations [rc»,:Bj] = 2i0^ n ' 



where f2,y is the generic element of the n-mode sym- 
plectic form 
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By definition, a quantum state p of n bosonic modes 
is said to be "Gaussian" when its phase-space Wigner 
function is Gaussian For this reason, a Gaussian 
state is fully characterized by its first and second-order 
statistical moments. These are the displacement vector 



Tr(xp) 



(3) 



2 



and the CM V, with generic element 

where {,} denotes the anticommutator. By definition, 
the CM is a 2n x 2n real and symmetric matrix. In order 
to be a quantum CM, it must also satisfy the uncertainty 
principle [20j 



V + ifl {n) > , 



(4) 



or an equivalent set of bona-fide conditions (for instance, 
see Ref. [2l[ for the case of two-mode CMs). In particu- 
lar, Eq. (j4j implies the positivity definiteness 



V > . 



(5) 



The simplest Gaussian state is the vacuum state, which 
corresponds to x = and V = I. 

Once that a state is prepared in a Gaussian states, 
its evolution can be such to preserve its Gaussian statis- 
tics. This is the case of Gaussian unitaries, which are 
defined as those unitaries transforming Gaussian states 
into Gaussian states. At the level of the second-order mo- 
ments, the action of a Gaussian unitary p — > U pU^ cor- 
responds to the congruence transformation V — > SVS T 
where S is a symplectic matrix, i.e., a matrix preserving 
the symplectic form SfJ^S 7, = fi("). A simple example 
is the beam splitter transformation. This is defined by 
the single-parameter symplectic matrix 



K(T): 



(6) 



where I is the 2x2 identity matrix and < T < 1 is the 
transmissivity of the beam splitter. In the Heisenberg 
picture, the beam splitter corresponds to the following 
Bogoliubov transformation of the quadrature operators 
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Finally, Gaussian measurements can be defined as 
those quantum measurements whose application to Gaus- 
sian states provides outcomes which are Gaussian dis- 
tributed [l| . When a Gaussian measurement is perfomed 
on a subset of modes of a bosonic system prepared in a 
Gaussian state, then the reduced state of the surviving 
(non-measured) modes is a Gaussian state. At the level 
of the second-order moments, the CM of the final state is 
connected to the CM of the initial state. As an example, 
when we homodyne one mode of a set of n bosonic modes 
in a Gaussian state, the formula of the final CM has a 
remarkably easy formula [l8l . [l9j . This is reviewed in the 
next section. 



III. COVARIANCE MATRICES UNDER 
HOMODYNE DETECTIONS 

Let us consider n + 1 bosonic modes in a Gaussian 
state. This (n + l)-mode Gaussian state p\ n has a CM 
that can be written in the blockform 



V„, = 



A C 
C T B 



(8) 



where A is the CM of the first n modes, B is the CM 
of the last mode, and C is a rectangular (2n x 2) real 
matrix accounting for the cross correlations. 



(n + 1) ■< 
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FIG. 1: An input Gaussian state pi n of n + 1 modes is ho- 
modyned in its last mode. The resulting output state p ou t of 
the first n modes is Gaussian. The input and output CMs are 
related by Eq. J9j . 

Now, let us homodyne the (n + l) th mode as shown 
in Fig. [I] performing the detection of the q quadrature. 
The output state p ou t of the remaining n modes is still 
Gaussian. In particular, this n-mode Gaussian state is 
described by the following CM [HGJ] 



a c(nBn) 



where 



n 



1 o 




(9) 



(10) 



If we detect the p quadrature, we have to consider Eq. © 
with the replacement 



n -> it 




1 



(11) 



Note that (IIBn) -1 is here a pseudoinverse (also known 
as Moore-Penrose inverse) since the matrix nBII is sin- 
gular. In general, for a rectangular and singular n\ x n 2 
matrix M, the pseudoinverse M _1 is an n 2 x n\ matrix 
G which minimizes the quantity 



:=£(H y ) 2 >0 



where H,*j are the entries of the matrix H := MG — I 
(with I being the n\ x n\ identity matrix). 

For the present problem, the pseudoinverse is easy to 
compute. In fact, if we set 



B 



bi b 3 



(12) 
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then we have 



and, therefore, 



nBn = 6, n 



(nBn) 1 = (6i) _1 n 



(13) 



(14) 



Similarly, for the detection of the other quadrature, we 
have 



so that 



n'Bn' = biri 



(n'Bn') 1 = (&0 n' . 



(15) 



(16) 



Since B > (being a reduced CM), then we have b\ > 
and the previous pseudoinverses of Eqs. (fT4)l and (|T6| are 
well-defined. 



IV. COVARIANCE MATRICES UNDER 
BELL-LIKE DETECTIONS 

Let us consider n + 2 bosonic modes in a Gaussian 
state. This (n + 2)-mode Gaussian state pi n has a CM 
that can be written in the blockform 



V- — 

* i.n. — 



C T B< J ) 



where A is the CM of the first n modes, 



B (2) 



Bi D 
D T B 2 



(17) 



(18) 



is the CM of the last two modes (labelled by 1 and 2), 
and 



C = (C X C 2 



c\ c k 2 



I k=l,n 



(19) 



is a rectangular (2n x 4) real matrix (describing the cor- 
relations between the first n modes and the last two 
modes). 

Here, we consider a Bell-like detection applied to the 
last two modes 1 and 2. This detection consists in apply- 
ing a beam splitter of transmissivity T, which transforms 
the input modes 1 and 2 into the output modes "+" and 
"— " , followed by two conjugate (q- and p-) homodyne 
detections as shown in Fig. [5] Thus, as a first step, let us 
apply the beam splitter symplectic matrix. The (n + 2)- 
mode Gaussian state p at the output of the beam splitter 
has CM 



V = [I (ri) ®K]V. m [I^©Kf J 



(20) 



(n + 2) ■{ 

P. 







►> «r> 



FIG. 2: An input Gaussian state pi n of n + 2 modes is subject 
to a Bell-like detection in its last two modes (labelled by 1 
and 2). The output state p out of the surviving n modes is 
Gaussian. The output CM is related to the input one by 
Eq. (T49)) . 



where K = K(T) is given in Eq. ([6]), and 

n 

l'")=0I (21) 



fc=i 



is the n-mode identity matrix (2n x 2n). This CM takes 
the blockform 



V = 



A C 

C T B( 2 ) 



where 



and 



C = CK i , 



§(2) = KB( 2 >K T . 



(22) 



(23) 



(24) 



More explicitly, the various blocks of the previous CM 
have the following expressions 



/ 



C = (Ci c 2 ) 



c{ c fc 



(25) 



' k=l,n 



with 



c 2 = -VI - TCi + Vrc 2 



and 



with 



B(2) = 



Bi D 

D T B 2 



(26) 
(27) 

(28) 



§! - TBi + (1 - T)B 2 + y/T(l-T)(D + D T ) , (29) 
B 2 = TB 2 + (1 - T)Bi - v/T(l-T)(D + D T ) , (30) 
D = y/T(l - T)(B 2 - Bi) + TD — (1 — T)D T . (31) 



4 



In terms of the previous blocks, the CM of Eq. (|22|) takes 
the more explicit form 



A Ci C 2 
V = | Cf Bi D 
Cl D T B 2 



(32) 



As already said, this CM describes the Gaussian state 
after the action of the beam splitter which transforms 
the last two modes 1 and 2 into the output modes "+" 
and 

Now, we apply the q-detection on the last mode "— ", 
and the p-detection on the next-to-last mode "+" . The 
detection of g_ implies the transformation of Eq. (0), 
which here reads 



A Ci 



V -> V' 

§)r(cf d t 



where 



r := (nB 2 n) 1 



(33) 



(34) 



In other words, the (n + l)-mode CM describing the first 
n modes and mode "+" is given by 



V 



a' a 
a T b' 



where 



A' = A c 2 r c 
B' = Bi D r D 



2 ) 
T 



and 



c = Ci - c 2 r d 2 



(35) 

(36) 
(37) 

(38) 



Now, let us apply the p-detection on mode "+". By 
using Eqs. (HJ) and ([TTj) . we get the final CM for the first 
n modes after the measurement, which is given by 



V -> v ouf = a' — c r' c' T 



where 



r' := (n'B'n') 1 



(39) 



(40) 



A. Simplification of the formula 

Here we simplify the formula for the final CM given 
in Eq. ([39]). Let us explicitly write the CM B^ 2 ^ of the 
detected modes by setting 



Bi 



fa fa 
fa fa 

D := 



, B 2 := 

^1 £3 
64 82 



fai fas 
fa 3 fa 2 



(41) 
(42) 



From these matrices, we can construct the symmetric 
matrix 



7 := 



7i 73 
73 72 



where 



7l := (1 - T)px + Tf3[ - 2y/T(l — T)5i 
12 := T/3 2 + (1 - T)0' 2 + 2y/T(l - T)S 2 



(43) 

(44) 
(45) 



and 



73 := VT(l - T)(fa 3 - /3 a ) - (1 
Then, it is easy to check that 

r-in, 

7i 



T)8 3 + T5i . (46) 



(47) 



and 



r' = 



71 

det7 



n' 



(48) 



Notice that previous Eqs. (|47|) and (|48|) are well-defined, 
since 71 > and det7 > 0, i.e., the matrix 7 is positive 
definite. Using Eqs. (|4"7) and (g5]), we can simplify the 
previous Eq. (|39[) . After some algebra, we get the main 
result of our paper, i.e., the input-output formula 



where 
Kn 



V nut = A - 



(1 



det7 



K 



22 — 



T"7 2 



v/T(l - T) 73 
T71 



-y/T(l-T) l3 

T _( - VT(1 - T) 72 

--T73 



-^(1-^73 
(1 - T)ti 

(1 



^)73 



v^a - t) 7 i 



(49) 

(50) 
(51) 

(52) 



Thus, the output CM V out of the surviving n modes 
after Bell-like detection is related to the input CM V in 
of the initial n + 2 modes of Eqs. (|17[) - (fiT)|) by means of 
the input-output relation of Eq. p9|) . where the matrices 
and 7 are completely characterized by the reduced 
CM B^ 2 ' of the detected modes and the transmission T 
used in the Bell-like detection. 

Note that, in Eq. (|4^]). the terms CiKyCj generate 
2n x 2n matrices, i.e., with the same dimensions of A. 
For instance, we have 




(53) 
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B. Specific cases 

Here we specify our formula of Eq. (|49|) to the case of 
standard Bell detection and that of heterodyne measure- 
ment. 



1. Standard Bell detection 

We achieve the formula for standard Bell detection by 
setting T = 1/2 (balanced beam splitter). In this case, 
we have 

li:=\{Pi+^)-h, (54) 

12 := \((32 + fc) + $2 , (55) 

and 

7 3 :=i(/33-/33-<5 3 +<54) • (56) 
Compactly, the 7-matrix takes the form 

7 = i(ZB 1 Z + B 2 -ZD-D T Z) , (57) 

where 

Z:=(J_° 1 ). (58) 
Also the K-matrices can be simplified. In fact, we get 

where 

Xi = ( a 1 ) , X 2 := f _ x 1 ) = n . (62) 

Since Ky = X^7Xj/2, we have that the formula of 
Eq. (|49|) becomes 

Vout = A " 2dcl^ E C,(Xf 7 X,)CJ , (63) 
where the 7-matrix is given by Eq. (|57l) . 

2. Heterodyne detection 

We achieve the heterodyne detection of the (n + l) th 
mode by setting T = 1/2 (balanced beam splitter) and 



considering the last mode in a vacuum state. As a matter 
of fact, this is equivalent to a standard Bell detection 
where the last mode is in a vacuum. Thus, we have 

B(2) = { B o 1) • C2 = ' ( 64 ) 

so that the global input CM is equal to 

v ™-( c A [B;) ffli - 

By setting B 2 = I and D = in Eq. ([ST)) , we derive the 
expression of the 7-matrix, which is given by 

7-i<nw + D-i( / U 1 flr *). W 

This matrix has determinant 

Q 

det 7 = -p 6»i := det Bi + TrBi + 1 . (67) 

Since C 2 = 0, the sum in Eq. (|6"3")l contains only the term 
with i = j = 1, i.e., 

V OMt - A _J_Ci(X? , 7X 1 )C?' . (68) 

Now, we can easily check that 

3^7X1 = ixf(ZBiZ + I)X! 

= - (n T Bin+l) = - ($7BiJ7 T +l) . (69) 
By using Eqs. ([57)1 and into Eq. (|55|) . we get 

V out - A id (nBi« r +l) Cf . (70) 

Similarly, if we heterodyne the (n + 2) th mode with the 
(n + l) th mode being the ancillary vacuum mode, we get 

Vout =A--j-C 2 (fiB 2 n T + I)C 2 " , (71) 

with 

0a := det B 2 + TrB 2 + l . (72) 

Note that this formula for the heterodyne detection co- 
incides with that reported in Ref. |1| (but where a proof 
was not explicitly given). 

V. CONCLUSION 

In conclusion, we have derived a simple formula for the 
transformation of CMs under Bell-like detections, where 
two modes of a bosonic system are subject to an arbitrary 
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beam splitter transformation, followed by homodyne de- 
tections. Specific cases correspond to the standard Bell 
measurement (where the beam splitter is balanced), and 
the heterodyne measurement (where the beam splitter is 
balanced and one input mode is the vacuum). Our for- 
mula can be applied in various contexts, including proto- 
cols of teleportation, entanglement swapping and quan- 
tum key distribution. In particular, it can be adopted 
to generalize these protocols to the case of experimental 



imperfections or asymmetric setups, which derive from 
the use of unbalanced beam splitters. 
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